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Abstract 



Starting from the gauge-invariant quark- ant iquark Green function the expres- 
sion for the nonrelativistic meson scattering amplitude on the external gluon field 
is obtained. The nonperturbative form of this amplitude is discussed. It is shown, 
that nonlocality of the nonperturbative vacuum field distribution plays an essen- 
tial role. As a concrete example, charmonium-nucleon scattering length is found 
and the possibility of the exotic charmonium-nucleus bound states is analyzed. 
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The structure and properties of the vacuum are among the most interesting 
questions of modern quantum theory. Of particular interest is the vacuum in 
nonabelian gauge theories. There are good grounds to believe, that this structure 
is very rich and complicated. Attempts to understand the hadron properties on 
the basis of the vacuum ones make an integral part of all papers in hadron physics. 

This paper is mainly devoted to the meson- meson scattering in QCD vacuum. 
In perturbative QCD quarks interact with each other by exchange of gluons and 
the role playing by one or other Feynman diagram depends on the kinematic 
region. It was also proposed a nonperturbative model for the scattering in the 
Regge limit s ^> m 2 , t ~ m? , where m 2 is a generic mass of the hadron (see papers 
0, P|) where quark- ( ant i) quark interaction is the result of QCD vacuum 
nontriviality and corresponding amplitude is proportional to the nonperturbative 
vacuum average of the gluon field stength 4-point function. 

Here and below we are going to describe the (nonperturbative) effects, which 
arise due to the interaction between nonrelativistic colour objects and nonpertur- 
bative gluon fields fluctuations in QCD vacuum in terms of the so called method 
of vacuum correlators (MVC) [Q, [15J. The corresponding physical picture in the 
case of nonrelativistic scattering is close to the one, developed in 0, ||, but 
there are some differences. First, for the ultrarelativistic Regge scattering the 
typical distances in the impact parameter plane are quite large and this is a rea- 
son to consider 'soft' Pomeron as a nonperturbative object. At the same time its 
structure may be more complicated, than it is proposed in ([]J - [[J), in partic- 
ular, it may include bound states, propagated in t-channel. Here we are going 
to discuss the role nonperturbative vacuum fields play in the heavy-quarkonium 
nonrelativistic scattering and to analyze the possibility of existence the exotic 
charmonium bound state in nuclei, which was suggested in |6[]. In this case we 
are interested in point-like type of interaction because of small quarkonium ra- 
dius and that's why we can apply the framework, which was discussed in H - 
. Second, it should be noticed, that we are looking for the scattering length for 
relatively slow heavy particle, which is propagated through the heavy nucleus. So 
we are going to explore usual nonrelativistic expansion over v/c, not the eikonal 
approximation, as it was done in ([[]]] - ||), and to restrict ourselves by the so 
called bilocal correlator, as one often does in the stochastic vacuum model [|J. 
The resulting expression for the amplitude corrects the formula, suggested in @ 
and provides the clear way for including the nonperturbative effects. 

A reasonable basis for the analysis may be, in our opinion, the hadron- 
hadron scattering description in terms of the connected Green function (G.f.) 
cluster expansion. As the first and the simplest step we are going to discuss 
meson G.f. only in this paper. 

As it was noticed, our technical basis is vacuum correlators formalism. In 
MVC vacuum is characterized by a set of cumulants like: 

(*) < F^(x 1 )$(x 1 ,x 2 )F p(j) (x2)...$(x n _ 1 ,x n )F^s(x n )®(x n ,x 1 ) > A 
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where F py (x) is the gluon field strength, $(xi,x 2 ) = Pexpig J A^dz^ is the 

parallel transporter and n = 1,2, ... . Also < ... >a means the average over all 
field configurations with the standard (euclidean) mesure. 

Consider the gauge-invariant Green function for a meson in an external colour 
field. It is easy to show that in bilocal approximation ||: 

= e- m ^~^ ■ GS(i2i(T0; R 2 {T 2 )) ■ G^Tx); f 2 (T 2 )) + 

+(H 2 f / dz, J dz' v < b«(z) b-(z') > A J d 3 R(r) j d 3 R(r>) 

° c c 

J d 3 r(r) J d 3 f(r')e- m( - T - T ^ ■ G^R^)- R(r)) ■ ^(^(Tx); f(r)) ■ 
. e -m(r'-r) . G\{R{t)] R(t')) ■ G^(r(r); f'(r')) ■ 
.e-miTz-r') . G^R(r');R 2 (T 2 )) ■ G$(f'(f); r 2 {T 2 )) (1) 
Here Tr{b^b v ) = | 6°6" and nonrelativistic G.f. : 



G^nin); f 2 (r 2 )) = i- J Dfe 



2/i 

r(ri)=ri 

The 6^ here are the external gluon fields, may be nonperturbative. The above 
Green function containes, up to the approximations have been made, all informa- 
tion about meson dynamics in the external field. In || this expression had been 
averaged over all vacuum fields and after it nonperturbative energy level shifts 
were extracted. We will understand these fields as belonging to the nonpertur- 
bative source, nucleon in our case, and the average over the gluon field states of 
this source will be made. The potentials V\ and V§ are the ones in the singlet 
and octet meson state respectively. We are going to analyze heavy quarkonium 
scattering, so for V\ and V 8 we limit ourselves by pure perturbative expressions: 

_ C F a s a s 1 _ 4 

Notice, a s is renormalized here on the characteristic quarkonia sizes - see discus- 
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It is convenient to express G^ in terms of fields operators. The simplest way 
to do it is by nonabelian Stokes theorem ||: 

J dz p J dz' u b^(z) b v (z') = J da^ p (w) J da uX (w') ■ 

c c 

■ ($(x , w) G„ p (w) x ) ■ ®(x , w') G vX {w') $(V , Xq)) 

Here da pu = a^drd/3 ; G pp {w) - gluon field tensor. The surface is parametrized 
by straight-line ansatze: w p {r) = z M (r)(l — (3) + z^{t)(3 and also 

/ dwu dw u dwn dw v \ 



d(3 dr dr d(3 
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Doing nonrelativistic expansion, it is easy to obtain: 

T 2 1 , 

J da w (w) G w (w) = J dr J dp I (Ef(r)) + (HL)- [p - ^) 
Ti o mm 

where = e^fc r i \ ~§P- an d E and i? are chromoelectric and chromomagnetic 
fields respectively, they are in the so called Feynman-Schwinger gauge here ( fixed 
point is x ), and in this gauge x ) = 1. 

Stricly speaking the last term in brackets is outside from our accuracy because 
it is of the order of ^. Deriving G^ we have omitted all such effects (the nonladder 
perturbative exchanges, spins of quarks et al.) Nevertheless, there may be some 
special cases, for example quarkonium-nucleus scattering via chromomagnetic 
fields, arizing from the nucleons Fermi motion, where nontrivial contribution 
from the H - term exists. 

G.f. in momentum space is as follows: 

G$(p in , Pout ) = J e-^dZRMe^^dZRidTz- 

■ < out\G2{z\ ) Z2)\in >, 

— * 

where Zi = (i?j(Tj), Tj), i — 1,2, and \in > ; \out > -are the asymptotic states of 
the singlet Hamiltonian. 

Also after returning to the M inkowskii space one can define the " amputated" 
G.f. (factor 2yU is a consequence of our nonrelativistic normalization ): 

G^iPin, Pout) = 2/x • lim (p] n - m 2 in ) G${p in , p out ) (p 2 0Ut - m 2 out ) 

Using spectral decomposition we obtain the expression for the nonrelativistic 
G.f. G^ in terms of E a operators: 



l l 



Gtf (Pi», Pin + '!) = 2m .(=£)j*frJ*frfdPSdiy 



< p\rE a (| + k) e _if (^- /J )(M|m > (2) 



Taking into account p in — (m + + E\ n , mv ), for nonrelativistic regime and 
also mi = m 2 = ^ one obtaines for q = : 



< aut\r'E\-k) e*"ti-r) \p > ko+E J_^ ■ 

■ < p\fE a {k) e-^ih-P) \i n > (3) 
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Up to a normalization factor G 2 is a scattering amplitude here. Formulae (2), 
(3) generalize a well-known result from 0,0. They have been obtained here in 
a gauge- invariant way, starting from the quark G.f. In an effort to demonstrate 
a nonabelian essence of (2), (3) it should be allowed the necessity to average the 
operator E a E a over the field source states. Let us define E a (k) = J d A xe %kx E a (x) 
and 

J d (^) e<^) < g 2 E?{x)E«{y) > A = 5 %] A q (x - y) (4) 

For zero momentum transfer: 

A q=0 (x-y) = A(x-y) 

The bilocal correlator A(x — y) = A(x — y, x 4 — y 4 ) is defined in an Euclidean 
space, hence one needs in Wick rotation k — > —ik^ ; it — > r: 



< 0Ut W - tk ; + Ei--H 8 U n \ m > ( 5 ) 

where A(x — y, k^) = J d(x± — y±) e %k ^ x ^-y^) A(x —y). Is it possible to consider 
the operator in brackets as a local one , as one often does? If A(r, k^) falls 
as a function of k^ on a distances about (T 5 ) _1 , characteristic values of k^ in 
the integral are controlled by [It = C^) -1 - Perturbative calculations give us 
E n = — ^ C 2°!^ , so for charmonium |e| ~ 600 Mev (e is a bound energy here, 
fi - reduced mass, for mi = m 2 = m/2 one has /i = m/4), a s is an effective 
strong coupling constant - with radiative corrections of the first order (see in 
H ) . Such distances are more than quarkonium size and more than T g , which 
are known from the lattice calculations and heavy quarkonium energy shifts is 
about (lGev) -1 . We will neglect x — y in A(x — y, fc 4 ), i.e. use approximation 
exp (ikf (1/2 — /3)) m 1. 

Nevertheless, even so 'space locality' exists, (with the rising accuracy if quarko- 
nium radius decreases) , there is a nontrivial dependence on (x 4 — y 4 ), in the 
correlator. A neglection of it is equivalent to the limit fix — > 0, which we shall 
discuss later on. So, it is obvious the necessity of special analysis of E(x)E{y) 
locality hypothesis both perturbative and nonperturbative cases. 

Let us make such analysis for charmonium-nucleon scattering with low energy. 
They usually use different parametrization for correlator A(x — y) ( 0, [[□]] ), 
taking into account lattice data or some models. We propose to use two forms 
for correlator A(x — y): 



A\{z) = A*(0) ; A h 2 (z) = A h 2 (0) e 



z 2 
To 2 



The letter "h" in 2 ( z ) reminds us, that the average here should be over all 
field configurations inside the hadron (not over the vacuum fields!), so there is a 
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question at this point, does the structure of these correlators coincide with the 
vacuum one? The positive answer to this question can be explained on the follow- 
ing way. Hadron parameters in the stochastic vacuum model are caused by the 
vacuum structure and the cluster expansion of the gauge-invariant Wilson loop 
depends only on the integral moments of the field correlators, so they contain all 
information about QCD vacuum, which may be investigated by colourless objects 
|3J . The way these moments are composed is not universal and it depends on the 
"hadron" (form of the loop, kinematics) but not the structure of correlators (see 



also |TT|). The only thing, that should be changed in our case is the normalization 



of the bilocal correlator. In the chiral limit (we use normalization condition from 



H) 



m N = < h\AQ^ \h > AO^ = 9 W - < 6 MM > A 
= 4^ CTL(x)CPL(x) (3(a s ) = - A a 2 s (6) 



where 6 = 9 for three flavours. From (4) and (6) it is easy to obtain the normal- 
ization condition for the correlator: 

Z \ 2 



A?, 2 (0)= i^-j 7r 2 -m N 

It is seen from the definition, that vacuum average of the correlators has been 
already subtracted here, so, we are indeed considering the scattering amplitude. 

In such a manner, amputated G.f. is proportional to an amplitude, as already 
noted. Normalized relativistic amplitude T is connected with G^ on a simple 
way T = 2m at • G^- Then we get: 

Ti = 2m N ■ — — / — < out\r/ 1 



2N C J 2tt 1 -ik 4 -Q + iO 

where operator Q = H$ + | E™ \ . Also 



A{(0,Jfc 4 ) = / dre ikiT A(0) e 

Finally, 



jii 

_ ~Tg 

OO 



T 1 = mm N 2 [-) n z < outlr 4 : , rAin > ; 

V3/ 1 H 8 — E\ n + T g — iO 1 

Analogously for Gaussian correlator A|(^) : 



T 2 = m niN 2 1 - ) 7r z < outWiTq^— I 1 — $ 



2\ 4 o I TJH%-E\ r - 



' 2 \~ ' \ 2 

T g 2 (W 8 -Bf") 2 



Ti\in > (7) 
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here = ^ /„* e~ l dt. 

For numerical analysis we choose m = M Vc = 2980 Mev, = 938 Mev. 
The limit T g — > oo (so called Leutwyler-Voloshin limit) leads to the regime the 



results of the paper |T2| was obtained in which. Indeed, it is easy to see, that 
using right value for the matrix element < 1 s | r,; H J_ El r i\^-S >, 



1 i 1 ii 7 i 

< Is n r» Is >= - ma 

1 H 8 — Et 41 



our results coincide: 



T = 7M ))c m c -(ym w aj (8) 

If one estimates the corresponding scattering length with the value a from []T2j 
a = 0.8 (Gef )~\ which reproduces the quarkonium bound energy for Coulomb 
wave functions, it is easy to obtain: 

I = — ; r = 0.012 fm 

8tt(M Vc + m N ) 



instead I = 0.05 fm in [T^. Taking into account the finite value of T g , namely 
(Tg) -1 = 1 Gev one has: 

h = — i — — r = 7 • 10" 3 fm 

87r(M„ c + m N ) 

T 

h = — ; r = 8 • 10~ 3 fm 

87r(M Vc + m N ) 

This result means, that charmonium-nucleus bound states are possible only for 
very heavy nuclei. Quarkonium-nucleon bound energy is so small, that potential 
wall deep enough may be only for A ~ 200. 

The above amplitudes equal to zero in another interesting case T g — > (mjv 
is fixed). In actuality, as it is evident from the foregoing, nonperturbative inter- 
action is of the different type in this case. 

As it is known, the system dynamics in this regime is natural to describe 
by the local potential . It means, it would serve no purpose to present the 
scattering in terms of gluon exchanges. In reality only bound colourless object can 
propagate over large distances in the impact parameter plane - it is reasonable 
to believe in this scheme, that it is a reggeized glueball. This possibility was 
supported sometimes ago by glueball Regge trajectory estimation |TJ]] in the 
framework of MVC, it seems close to the phenomenological one for the Pomeron. 
The next step should be the Regge scattering description in terms of QCD vacuum 
values, taking into account also direct nonperturbative interaction according to 
and perturbative exchanges in the nonperturbative vacuum, where they are 
important. This work is in progress now. 
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